function isCorrectDirection = get_reaction_dirrection_by_thermodinamics(RKEGGID,equation, CKEGGID)
    thermodinamicsValues = evalin('base', 'thermodinamicsValues'); 
    defines = evalin('base', 'defines'); 

    isCorrectDirection = false;
    index = strmatch(RKEGGID, thermodinamicsValues);
    isOnLeft = false;
    value='0';
    if ~isempty(index)
        if ~isempty(thermodinamicsValues{index,2})
             value = thermodinamicsValues{index,2}(1,:);
        else
            isCorrectDirection = true;
        end
    else    
        isCorrectDirection = true;
    end
    if ~isempty(value)
        value = str2num(value);
        if abs(value) > str2num(defines.thermodinamicTreshold)
            if regexp(equation(1: regexp(equation, '<=>') - 1),CKEGGID) > 0
                 isOnLeft = true;
            end
            if strcmp(defines.isSubstrateDirection,'true')
                if (value < 0 & ~isOnLeft) | (value > 0 & isOnLeft)
                    isCorrectDirection = true;
                end
            else
                if (value < 0 & isOnLeft) | (value > 0 & ~isOnLeft)
                    isCorrectDirection = true;
                end
            end
        else
            isCorrectDirection = true;
        end
    end
    clear defines;
    clear thermodinamicsValues;
end